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Abstract 

Linear conductance of junctions formed by graphene flakes with order of nanometer-thick 
electrodes attached at the corners of the flakes is studied. The explored structures have sizes 
up to 20 000 atoms and the conductance is studied as a function of applied gate voltage varied 
around the Fermi level. The hnding, obtained computationally, is that junctions formed 
by armchair-edge flakes with the electrodes connected at the acute-angle corners block the 
electron transport while only junctions with such electrodes at the obtuse-angle corners tend 
to provide the high electrical conductance typical for metallic GNRs. The finding in case 
of zig-zag edges is similar with an exception of a relatively narrow gate voltage interval in 
which each studied junction is highly conductive as mediated by the edge states. The contrast 
between the conductive and insulating setups is typically several orders of magnitude in terms 
of ratio of their conductances. Main results of the paper remain to a large extent valid also 
in the presence of edge disorder. 
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1 Introduction 


Extensive research of graphene was initiated by its nniqne electronic properties and at the same time 
by the possibility to prepare graphene samples in laboratories [1, 2], On the theoretical side, the 
dc electronic transport in ideal graphene and GNRs is the process of fnndamental interest. Experi¬ 
ments addressing electronic transport in graphene necessarily employ hnite graphene flakes, nsnally 
rectangnlar graphene nano-ribbons (GNRs). Therefore, in GNRs, their hnite spatial dimensions 
and conseqnently the bonndary conditions mnst be taken into acconnt. It is known that armchair 
GNRs (AGNRs) can either be metallic or semi-condnctive, depending on their width [3, 4]. The 
zig-zag GNRs (ZGNRs) are always metallic and snpport special localised states at their edges [3]. 
Gontacts of the electrodes to a GNR are often assnmed to extend across the entire width of the 
ribbon. While the condnctance properties of snch jnnctions have been freqnently stndied [5, 6, 7] 
(see also review [2] and references therein), little is known abont graphene hakes contacted at their 
corners. In their compntational work [8] Gnong et al. reported the absence of edge states near 
the 120° corners of ZGNRs. This hnding is conhrmed also by onr calcnlations and will be refer¬ 
enced below in the text. In another compntational stndy Bornnda et al. investigated cnrrent how 
throngh rectangnlar GNRs assnming electrodes contacted at certain positions at opposite edges of 
the GNR [9] . Based on their simnlations the anthors were able to diherentiate ballistic and dihnsive 
regimes of transport in graphene. Related to onr present work is also the theoretical analysis of 
a dnal-probe STM setnp on graphene by Settnes et al. [10]. The anthors compnted condnctance 
maps assnming the two STM tips at varions positions on graphene snrface, inclnding non-ideal 
samples. We however conld not hnd any analysis of graphene jnnctions with order of nanometer 
thick electrodes attached at the corners of the hake. Therefore the goal of this work is to present 
onr resnlts on the dc linear condnctance of graphene hakes of varions shapes and edges and with the 
electrodes connected at the corners of the hakes. Particnlar investigated samples have the shapes 
of trapezoids, triangles, rhombi and rhomboids, all with angles at corners being integer mnltiplies 
of 30°. We hnd that the condnctance strongly depends on the angles at the corners to which 
the electrodes are connected. Acnte-angle corners (30° and 60° stndied here) in most cases do 
not enable good electronic transport. The opposite is trne for electrodes attached at obtnse-angle 
corners (120° stndied here): they often provide high condnctance levels typical for metallic GNRs. 
We describe GNRs’ electronic strnctnre using a tight-binding (TB) approximation with hoppings 
up to third nearest neighbor (NN) as specihed below in section 2. The model of the electrodes 
uses the hrst NN approximation and will further be described as well. Main results of the paper 
are explained within section 3, considering graphene flakes with shapes of symmetrical trapezoids 
and triangles. Section 4 supports the main results by a local density of states analysis. Section 5 
provides a summary and discussion of the results. Gomparisons to results obtained within the hrst 
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Figure 1: Schematic graphical representation of our model of the electrodes and their coupling to 
the graphene flake. In this example the flake has zig-zag edges and the electrode is attached to its 
60° corner. The electrode itself is formed by a bunch of monoatomically thin semi-inhnite mutually 
decoupled wires. In typical setups we use 35-36 wires per electrode. Here we show only 12 of them 
for easier visualisation. The interatomic distances within each wire are a, i.e. the same as is the 
NN distance in graphene. See also appendix B in which we consider a model with wires mutually 
coupled along a hnite length. 


NN TB approximation (INNTBA), further details and examination of the model of the electrodes, 
and a section describing results for samples with edge disorder are supplied as appendices. Supple¬ 
mentary information (SI) contains results for graphene flakes of several additional shapes including 
a system with wide contacts to electrodes, further results within the INNTBA and a section with 
data for the flakes terminated by imperfect edges. 


2 Models and Methods 

Our model of GNR’s electronic structure employs one explicit electron per atom, the independent- 
electron approximation and a TB Hamiltonian with hoppings up to the third NN included (3NNTBA) 
According to [11], we set the hopping parameters to values t-Q = —2.97eV, tg = —0.073eV and 
fg = —0.33eV. Magnitude of te is conveniently used as the energy unit. Interatomic overlaps are 
neglected. The edges of the flakes are assumed to be terminated by hydrogen atoms electrons of 
which do not explicitly enter the TB Hamiltonian; see for instance [12]. Although the TB mod¬ 
elling is a relatively basic-level one, it has, even in the INNTBA, in many cases been proven to 
provide satisfactory description of conductance properties of graphene structures on at least semi- 
quantitative level which is sufficient for the present study. Such models are frequently employed in 
theoretical studies of graphene samples [7, 13, 14] and carbon nanotubes [15]. 

In the present study the whole system including the electrodes is composed of identical atoms. 


3 











Each of the two electrodes in our model is composed of a number (from tens to hundreds) of 
monoatomically thin mutually non-interacting wires (hgure 1). Within each wire the INNTBA 
with the parameter t-^ is used to describe the electronic structure. This simple model of the 
electrodes was found very convenient in [16] where it was used for GNRs and provided results in a 
semi-quantitative agreement with different models of the electrodes [17, 14], The coupling of the 
wires to the GNR is again described by the hrst NN model using the parameter t-Q. See appendix B 
where we provide an examination and justihcation of such model of the electrodes for the use in 
the present work. 

Each atom of the system is characterised by its on-site energy. In equilibrium and in absence of 
any external held these energies are set to zeros for all the atoms of the entire system. The effect of 
the applied bias voltage U as well as of a gate voltage are modelled by variations of the on-site 
energies of the atoms in the electrodes and/or in the graphene hake [16]. Variation of the gate 
voltage represents the shift of the chemical potential in the graphene hake away from the neutrality 
point and is considered within a 0.2 I^bI wide interval around the Fermi level. We apply almost 
limitingly small (still numerically hnite) values of the bias voltages t/, compute resulting dc currents 
I through the junctions and consequently obtain the linear conductances Gun = /U. The 

considered model can now be briehy summarised: Hamiltonian of the entire system is 


H — ^ Hi^i! a\aii , 
1 , 1 ' 


( 1 ) 


with I and I' running over atomic sites of the entire system (including the electrodes) and a|, 
a;/ being the creation and annihilation operators of an electron in the TB orbitals |/) and |/'), 
respectively. It is assumed that the TB orbitals form and orthonormal basis set: (/|/') = 5^/. The 
matrix elements of the Hamiltonian in our model are set to 


Hi,v = 


eU 

1 = 1' within the wires of the left electrode 

0 

1 = 1' within the wires of the right electrode 

(eVg + eU/2) 

1 = 1' within the flake 

tB 

/, V being any hrst NN 


/, V being second NN within the hake 

t'B 

Z, 1' being third NN within the hake 

0 

all other 1, V 


( 2 ) 


Since the monoatomically thin wires are assumed to be mutually decoupled, matrix elements Hi^i/ 
with I, V belonging to diherent wires are all zeros, even if the two wires belong to same electrode. 
The exceptions are the examples used in appendix B in which we examine a more complex model 
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of the electrodes. As mentioned above, we explore the linear regime, choosing particular numerical 
value of the bias voltage U = 10“® Itnl/e. Around and below this magnitude of U the I{U) function 
is perfectly linear; even higher values would still provide the linear regime. The contribution ef//2, 
although applied within the graphene flake [line 3 in equation (2)], does not have a real effect in 
this particular study because of the small value of U. We note in passing that application of the 
INNTBA in the entire system would yield a perfectly symmetric functional dependence of Gun: 
Giin(— 14 ) = Giin( 14 ), thanks to the symmetry of the dispersion relations of the wires as well as 
of graphene [1]. Inclusion of the hopping matrix elements up to the third NN in graphene flakes 
breaks the symmetry. 

In addition, the extension of the TB hopping range in the graphene flakes causes, among other 
features, a horizontal shift of the calculated Giin(I4) prohle compared to the results from the 
INNTBA. This shift results from the approximate model (2). Its value would change with in¬ 
clusion of overlaps in the TB model. To facilitate presentation of the results we apply a proper 
compensating horizontal shift AI 4 = +0.072 Itnl/e to all plots obtained by the 3NNTBA so that 
they match the symmetry of corresponding INNTBA curves as much as possible in cases of the 
armchair-edged structures. We apply the same compensating shift AVg also for zig-zag-edged 
flakes although in this case the plots can not be made even approximately symmetric. Hence our 
horizontal axes description in terms of quantities formally equals to I 4 + AI 4 . 

Electron currents in the present work are obtained using a Landauer-type formula in the same 
way as we used in stationary calculations in [16]. Because an electrode in our work is modelled as 
a bunch of mutually non-interacting monoatomic wires, we formally have a multi-terminal system 
and use a multi-terminal formalism. For example, each of the junctions depicted in hgure 2 uses 
two 36-wires thick electrodes. Such junctions are formally treated as 72-terminal systems within 
our approach. We provide more details on the model of the electrodes in appendix B where we 
also demonstrate that the model based on the non-interacting wires does not yield any signihcant 
consequences to the calculated linear conductances. The core part of the multi-terminal formalism 
are the quantum-mechanical calculations which employ the scattering approach with wavefunctions 
obtained from Green function formalism [18], again in the same way as we used in [16]. A scattering 
approach and the Landauer-type formula were used also in our work [19] for systems with just two 
terminals; see equations (24) and (25) therein.^ Here we use the multi-terminal generalisation of 
the method and instead of the simple scattering approach we employ the more general one based 
on the Green function technique which allows us to obtain the wavefunctions. 

^Formula (24) in our work [19] has the density-of-states factors included in the definitions of coefficients A and 
C. 
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3 Results 


Conductance properties of a graphene-flake formed nano-junction depend on a number of parame¬ 
ters and conditions: on the intrinsic properties of the flake itself (its shape, size) and also on the way 
how and where the electrodes are attached to the flake. Among the plethora of possible scenarios 
we focus on several representative ones with the aim to study junctions with the electrodes at the 
corners. We present our results in hgures which display particular considered structures together 
with the linear conductance Gun as a function of the chemical potential variation (here represented 
by the gate voltage) around the neutrality point. Results for a limited set of structures only are 
described in the text, which however should be sufficient to explain our main hndings. Results for 
other structures can be found in the SI. In our description of the physical structure sizes we often 
employ the graphene NN distance a ~ 1.42 nm as well as its lattice parameter b = ay/3 ~ 2.46 nm. 

3.1 Armchair edges 

As representative structures with armchair (ac) edges we consider symmetrical trapezoids with 
acute angles of 60° at their base (left panel of hgure 2). These flakes support four different basic 
attachments of the electrode pair: (i) both electrodes at the acute angles (AA-AA setups), (ii) one 
electrode at the acute angle and the other at the nearby obtuse angle (AA-OA setup), (hi) one 
electrode at the acute angle and the other at the opposite obtuse angle (the other AA-OA setup), 
(iv) both electrodes at the obtuse angles (OA-OA setup). The top three structures in hgure 2 
share the same base length L = 176 a ~ 25.0 nm. The equilateral triangle in this context can be 
considered as a limiting-case trapezoid. The two trapezoids at the bottom represent wider and 
overall larger structures {L = 353 a ~ 50.1 nm) and allow us to observe size-dependence of the 
computed results. The heights of the samples are 14 b, 32 b, 88 b, 23 b and 46 b from the upper-most 
structure to the one at the bottom. Total numbers of atoms composing given hake are typed on 
the images. Each electrode considered for hgure 2 is 36 monoatomic wires thick. While possible, 
we choose the shape of the contact areas the same for corresponding corners.^ All trapezoids 
in hgure 2 have (vertical) widths corresponding to metallic AGNRs [3, 4]. Right panel of the 
hgure displays our results expressed in terms of the linear conductance as a function of the gate 
voltage -|- AV^. Importantly, we hnd most of the setups with both electrodes at the 60° corners 
(the AA-AA junctions) to be relatively insulating near the neutrality point (black solid curves in 
hgure 2) with their linear conductances Gun being several orders of magnitude below those for OA- 
OA setups (blue solid curves). The current-blocking behaviour is mostly retained also at elevated 
chemical potentials apart from the isolated peaks. On the contrary, conhgurations 120° — 120° 
(the blue plots) are highly conductive, with Gun ~ 0.2 — 0.3e^/h roughly independent on the gate 

^The largest trapezoid has sharper corners at its obtuse angles, compared to the smaller structures. 
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Figure 2: Results for the trapezoidal flakes with armchair edges. Left panel: graphical represen¬ 
tation of the flakes. Red-colored atoms in the corners are directly coupled to the electrodes. Each 
individual considered setup employs just two of the four colored corners serving as the electrodes- 
attachment areas. The electrodes are 36-monoatomic wires thick. The numbers of atoms in each 
flake are displayed by labels 3168, ..., 19 158. a ~ 1.42 nm is the NN distance in graphene. See 
main text for more details. Right panel: linear conductances Gun plotted as functions of the gate 
voltage for the structures shown in the left panel in the same vertical order. t-Q is the nearest- 
neighbor tight-binding parameter used for the entire system and e is the unit charge. Solid black 
curves: both electrodes at the 60° corners. Dashed red curves: one electrode at the 60° corner, 
the other at the nearby 120° corner. Solid green curves: one electrode at the 60° corner, the other 
at the opposite 120° corner. Solid blue curves: both electrodes at the 120° corners. Legends like 
120° — 120° apply globally within the whole figure. For graph (e) the two plots only have been 
calculated. The independent variable on the horizontal axes is V^-|-AV^, with AV^ = -1-0.072 Ifsl/e 
as explained in section 2. 


voltage, with the exception of the largest flake [figure 2(e)]; even the latter example provides a huge 
contrast (4-5 orders in magnitude) between the conductances of the OA-OA and AA-AA setups at 
the Fermi level. The order of 10“^ e^/h conductance magnitudes are typical for metallic graphene 
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nano-ribbons [6, 7]. The case of the mixed setups (one electrode at the 60° angle while the other 
at one of the 120° angles) is more complicated; see dashed red and solid green plots in figure 2. 
Typically, linear conductances of these AA-OA setups take intermediate values. An interesting 
feature in some of the AA-OA setups, found for example in figure 2(b), is that the prohle of the 
Giin function is mostly independent on the choice of the particular pair of the AA-OA attachment 
corners: both the nearby-corners case and the opposite-corners case yield similar curves; compare 
the dashed red and solid green plots on the figures. The similarity is reduced with increasing 
relative length of trapezoids. This is an intuitively comprehensible feature because for long-narrow 
trapezoids the mutual positions of the contacts in the two 60° — 120° setups are quite different. In 
SI we show that there is no such similarity in case of ac-edge terminated long-narrow rhomboids. 
Note that the AA-OA junctions have not been calculated for the largest trapezoid, hgure 2(e). 
Finally we refer to appendix section A which provides results employing the linear scale on the 
conductance axis [figures 8 and 9]. 

3.2 Zig-zag edges 

For graphene flakes with ZZ edges we perform an analysis analogous to that in previous section, 
choosing now electrodes 35 monoatomic wires thick. We again choose structures of trapezoidal 
shapes. The samples displayed in hgure 3 are of two different lengths, 100 b and 204 b. Their heights 
from the upper-most hake to the lowest one are 24.5 a, 56 a, 150.5 a, 41a and 80 a. Generally valid 
diherences if compared to the ac-terminated hakes (hgure 2) are that (i) The AA-AA junctions 
now provide a signihcant conductance within a 0.025 Itel narrow gate voltage window; this is an 
ehect of the special localised ZZ-edge state and can be found signihcant if TB hoppings up to 
the third NN are included; this high-conductance regime would not be found for AA-AA setups 
within just the INNTBA as we have checked and provide several comparisons in appendix A, most 
relevantly in hgure 11(b). (ii) The conductance ratios of the OA-OA setups and the remaining ones 
are now larger for most of the studied Vg range, (iii) The Gun spectra are signihcantly asymmetric 
around the Fermi level, (iv) Rapid oscillations are visible in the Giin(I4) spectra within the high 
conductance window. 

As mentioned above, the high conductance regime of the AA-AA setup found in the vicinity 
of the Fermi level comes from the contribution of the ZZ edge state, hgure 3 shows the rapid 
oscillations of Gun in this gate voltage range. The oscillations can be more clearly seen in data 
obtained within the INNTBA and consequently more easily interpreted. We demonstrate such 
results in SI for long-narrow rhomboids. The presence of the oscillations even at gate voltages 
arbitrarily close to the zero can be understood taking into account the ehect of the electrodes 
and coupling of their modes to the hake: the semi-inhnite electrodes support modes at energies 
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Figure 3: Results for the flakes with zig-zag edges. All description as well as the color coding is 
analogous to that in hgure 2. The electrodes are 35-monoatomic wires thick now. Here as well as 
in all other graphs displaying the conductances obtained within the 3NNTBA model we apply a 
uniform horizontal shift AV^ = -1-0.072 Itnl/e to all computed plots as described in section 2. The 
distance b is the lattice parameter of graphene given by the relation b = ay/3 ~ 2.46 nm, where 
a ~ 1.42 nm is the NN distance in graphene. 


arbitrarily close to the Fermi level. 

Away from the energies or in absence of the ZZ-edge induced peaks, the electronic transport 
between the corners is predominantly mediated by the bulk modes, not by the edge states. This 
is demonstrated by the 120° — 120° plots (the OA-OA setups) in hgure 3 which shows the almost 
uniform high conductances across the whole interval of gate voltages. In this way our data in 
terms of the conductance rehect the fact that the ZZ edge state does not exist at and close to 120° 
corners [8]. The transport between the 120° corners of the ZZ-edge structures is only marginally 
affected by the ZZ edge states. 

At the end of this section we conclude that outside the narrow range of the ZZ-edge conductance 
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Figure 4: The linear conductance for a junction formed by an isosceles-triangle shaped graphene 
flake with the electrodes attached at its corners. The inset shows an image of the flake. The size of 
the structure is 10 094 atoms. The base edge is of the ac type while the edges along the sides are of 
the ZZ pattern. The electrodes’ thickness is 33 monoatomic wires. Other description is analogous 
to that in previous two figures. 

window the linear conductance of ZZ-terminated trapezoid is typically vanishing unless the two 
electrodes are attached at the OA corners. Outside the ZZ-edge conductance window the contrast 
between the conductances of the OA-OA and AA-AA setups is several orders of magnitude. The 
other similarity to the case of ac-edge terminated trapezoids is that the results for the two AA- 
OA setups are very close each other (red and green plots in figure 3). Contrary to the ac edges, 
the proximity of the two AA-OA curves is now found also for long-narrow structures; additional 
examples of this feature can be found in SI on the case of rhomboidal flakes. Although we do 
not analyse this effect we can say that it is related to the electronic modes which are used for the 
transport: the edge states and the bulk states. Finally, we again refer to appendix section A, in 
particular hgures 10 and 11, which use the linear scale on their conductance axes. 

3.3 Corners formed by different pair of edges: isosceles triangles 

Besides the trapezoids, we have addressed also graphene flakes of several other shapes, including 
rhomboids, rhombi and equilateral triangles. All of these feature 60° and 120° angles at their 
corners. Not surprisingly, they provide basically the same picture as the results for trapezoids 
described above. Some of these results are show in SI. As a different example here we consider a 
flake of an isosceles triangle shape with 30° acute angles, the ac edge at the base and the ZZ edges 
along its sides. In other words, in the 30° angle corners the two crossing edges are of the different 
types (ac and ZZ). Based on the above hndings, we intuitively expect a large contrast between the 
conductances of the setups 30° — 30° and 30° — 120°. Quantitative results shown in hgure 4 conhrm 
this expectation, with the exception of the isolated resonances and, more significantly, with the 
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Figure 5: The linear conductances for different sizes of the contact areas (marked by red color). 
The equilateral triangle-shaped flakes have ac edges. Each triangle is of the same size, consisting 
of 10 620 atoms. The numbers 36, 126 and 468 (both in graphs and at the triangles) denote the 
number of contacted atoms per corner, i.e. the electrode’s thickness in terms of the number of the 
composing wires. 

exception of the 0.025 {tsl narrow high-conductance channel which is the signature of the presence 
of the ZZ edge, although now only single one at each acute-angle corner. We can expect that results 
for the case of a ZZ base would be similar to those in hgure 4. 

3.4 Contact size effect 

Naturally, the conductance depends on the size of the contact area which in our model is represented 
by the electrode’s thickness and is quantihed by the number of the monoatomic wires per electrode. 
We have seen that graphene flake junctions can be highly conductive especially for the OA-OA 
contacts even if the contact area corresponds to just 35 or 36 atoms. For these contact sizes the 
AA-AA setups have usually been found relatively insulating for most of the applied gate voltages, 
the exception being ZZ-edge mediated transport window. In this section we study if and how the 
predominantly insulating behaviour is modihed due to increased contact areas. We opt to study 
this effect on equilateral-triangle shaped flakes. The results are shown in hgures 5 and 6. For the 
ac-terminated structures we hnd that the increase of the contact areas leads to the broader peaks 
of the conductance spectrum while the peaks’ maxima remain roughly unchanged. The minima of 
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Figure 6: Same kind of results as in figure 5 but now for the flakes with ZZ edges. The triangle is 
composed of 10 401 atoms. The numbers 35, 143 and 360 denote the number of contacted atoms 
per corner. 

the spectrum become systematically higher for thicker electrodes. Although the conductance for 
the largest contact area (468 atoms) is signihcant at a whole range of elevated gate voltages, the 
ratio of the conductance to the number of contacted atoms remains small compared to the cases 
when contacts are made at the obtuse angles (hgure 2, the blue curves). 

In case of the ZZ-edge terminated triangles (hgure 6) the conductance displays the central 
ZZ-edge mode contribution and few additional very narrow resonances. The spectrum within the 
high-conductance window does not change signihcantly with the thickness of the electrodes. The 
remaining part of the calculated spectra in hgure 6, while being rather asymmetric in comparison 
to the ac case (hgure 5), exhibits a similar dependence on the electrodes’ thickness. 



4 LDOS analysis 

The current-blocking vs. conductive behaviour reported above can be understood by studying the 
electronic modes of the graphene hakes. We use the local density of states (LDOS) of isolated hakes 
as the tool. We start from the dehnition of the density of states of an isolated hake projected on 
an atomic site 1: 

P,{E) = J2\^l\i,i)p6{E-Et), (3) 
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Figure 7: Local density of states (LDOS) in corners of isolated graphene structures. The function 
is computed over the red-marked corner regions of two trapezoidal graphene flakes, (a) LDOSs 
in the 60° and 120° corners of the armchair flake that is shown on hgure 2, the second structure 
from the top. (b) LDOSs in the corners of the zig-zag flake that is shown on hgure 3, the second 
structure from the top. The legends referring to the angles apply globally within both graphs. The 
LDOS curves are plotted as they come from calculations, without any additional shift like AV^. 


with \l) being the atomic orbital at the site, tj^i is eigenfunction of the isolated hake and Ei the 
associated eigenenergy. We acquire LDOS information from a region of interest by summing up 
over corresponding sites: 

LDOS(E) = E f,(E). (4) 

/Eregion 

In the numerical evaluation the ^-functions are replaced by normalised gaussians 


1 


w 




exp 


(e-ea 

\^w) ’ 


(5) 


with w = 0.015 |tB|- The spectra are plotted in hgure 7. Its graph (a) shows LDOSs for the case of 
the ac-terminated trapezoid composed of 6324 atoms and displayed in hgure 2, the second structure 
from the top. Graph (b) of hgure 7 shows LDOSs for the ZZ-terminated trapezoid composed of 6306 
atoms and displayed in hgure 3, the second structure from the top. The /-summation in equation (4) 
runs over the red-marked atoms of the chosen corner [36 atoms in case (a) and 35 atoms in case 
(b)]. In the ac case the LDOS at the 120° corner (solid red line) is clearly larger than the LDOS 
at the 60° corner (dashed black line) within the central range of energies relevant for the electronic 
transport studied here. In the ZZ case [hgure 7(b)] the situation is more complicated because of 
the presence of the special localised edge state [3] at zero energy and the corresponding peak in 
the LDOS. Still, apart form the edge state contribution, the values of LDOS at the 120° corner 
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are again by far larger than the LDOS at the 60° corner. Contribntion of the edge state to the 
electronic transport between corners of the flakes studied here is specihcally limited or suppressed 
in some cases: this state does not exist at and close the 120° angles of graphene flakes as it was 
shown in [8] and as our analysis conhrms, see discussion in section 3.2. Our LDOS data for the ZZ 
case conhrm the absence of the edge state: the solid red plot in hgure 7(b) (corresponding to the 
120° angle) reaches a maximum of only about 5 on the given interval while the dashed black curve 
(the 60° angle) has its maximum at about 68 (beyond the axis scale). 

Finally we note that while in the graphs with Gun we have applied a uniform horizontal shift 
AI 4 = +0.072 Ifnl/e to calculated data (see section 2 for explanation), no such a shift is used for 
the LDOS data. 


5 Discussion and Conclusions 

In this work we studied stationary zero-temperature quantum transport through junctions formed 
by non-rectangular graphene flakes with about nanometer or several nanometers thick electrodes 
attached at the corners of the flakes. Such structures are perfectly compatible with hexagonal 
graphene lattice. Typical shapes are rhomboids, rhombi, trapezoids and triangles. They feature 
corners with angles being integral multiples of 30°. Many of these structures have all their edges 
of the same kind: either armchair (ac) or zig-zag (ZZ). Isosceles triangles with 30° acute angles 
provide combinations of the ac and ZZ edges. Studied flakes have sizes up to almost 20 000 carbon 
atoms. Assuming a low bias voltage regime we computed the linear conductance as a function of 
a gate voltage applied to the flake. The computational methodology was based on a tight-binding 
model with hoppings up to the third nearest neighbor included and on the scattering approach. 
The later technique was employed in stationary calculations in our recent work [16]. 

Main hnding of the present work is that the conductance strongly depends on the size of the 
angles at the corners to which the pair of the electrodes is attached. Especially for flakes with ac 
edges we can say that junctions with both electrodes connected at obtuse-angle (OA) corners provide 
1-5 orders of magnitude higher conductance compared to setups with the electrodes attached at 
acute-angle (AA) corners. The OA-OA junctions often provide high conductance values typical 
for metallic graphene nano-ribbons. On the contrary, the pair of the electrodes attached at acute 
angles forms a relatively insulating junction. The contrast between the conductances of OA-OA 
and AA-AA setups extends beyond one order of magnitude within a 0.15 \t^\ wide interval of the 
chemical potential tuned around the Fermi level, being the nearest-neighbor hopping parameter 
of the tight-binding model. Exceptions are a few narrow isolated resonances at which also AA-AA 
setups yield a high conductance. Quantitative values depend on the dimensions of particular flake. 

For graphene flakes with ZZ edges the situation is complicated by the existence of the well- 
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know edge state at the Fermi level. This mode couples with the electrodes and provides the high 
conductance also for AA-AA setups, but only within a 0.025 Itel narrow interval of the chemical 
potential close to the Fermi level. For the remaining part of the spectra in the ZZ cases we again 
hud a huge contrast between the conductances which can be even 1-2 orders of magnitude larger 
than for ac edges. The spectral range in which the high contrast is found is also significantly wider 
that in cases of flakes with the ac edges. 

We have found that the reported behaviour does not change qualitatively within a range of 
electrodes’ thicknesses, from sub-nanometer values to several nanometers. Another important fea¬ 
ture of real samples is the presence of edge defects and their impact on the electronic transport. 
It is known that especially AGNRs are more affected by the presence of the edge disorder. On the 
contrary, flakes with ZZ edges are much more robust in this respect thanks to the presence of the 
edge state [12], see also [20, 21]. Although we did not address this problem in a systematic manner, 
we have at least examined two representative trapezoidal flakes with imperfect edges, an ac-edge 
terminated trapezoid with random edge defects as well as a ZZ analogue. We provide details in 
appendix C and additional material in SI. Our findings fully confirm the aforesaid properties: espe¬ 
cially for graphene flakes with ZZ edges the effects reported in the present paper remain preserved 
to a large extent even in the case of signihcant edge disorder. In the case of ac edges we observe 
significant drop in the conductance of the junction which would otherwise be highly conductive with 
perfect edges. Still, even in the ac case, the conductance of the OA-OA setup remains to be several 
orders of magnitude larger compared to the AA-AA setup within a range of gate voltages. We 
expect that experimental confirmation of the effects reported here could be realised using proper 
carbon nano-ribbon or nanotube electrodes, possibly providing also a mechanical support for the 
graphene flake. 


Acknowledgements 

This work was supported in parts by the Slovak Research and Development Agency under the 
contract No. APVV-0108-11 and by the Slovak Grant Agency for Science (VEGA) through grant 
No. 1/0372/13. The structures were visualised using VMD [22]. The author thanks Peter Bokes 
and Peter Markos for stimulating discussions. 

A vs. 3*^^ nearest neighbor model 

All results in main text were obtained within the 3NNTBA. In addition we have performed several 
calculations in which TB hoppings up to only the hrst NN were included (INNTBA). In this section 
we compare results by INNTBA and 3NNTBA approaches and in this way evaluate the effect of the 
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Figure 8: Linear conductances compared for the INNTBA (solid green line) and 3NNTBA 
(dashed red line) models. The examined structure is the ac-terminated trapezoid composed of 
6324 atoms (hgure 2). The functions for the 3^'^ NN model have been horizontally shifted by 
AV^ = +0.072 Itnl/e as specihed in main text, section 2. Graph (a) shows results for the electrodes 
attached at the 120° angles. Graph (b) shows analogous results for the 60° angles. The models 
of the electrodes are identical to those used for results in hgure 2. See also hgure 9 for a detailed 
view on the low-V^ interval of the present graph. Note that for hgures within this section we use 
the linear scales also on vertical axes. 

extended hopping range on the conductance. The wires forming the electrodes are always modelled 
within the INNTBA. Examined systems are the ac-terminated trapezoid composed of 6324 atoms 


trapezoid, armchair edges, 6324 atoms, detail 
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Figure 9: A detailed view of the results from hgure 8. 
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trapezoid, zig-zag edges, 6306 atoms, overall view 
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Figure 10: Linear conductances compared for the INNTBA (solid green line) and 3NNTBA 
(dashed red line) models. The examined structure is the ZZ-terminated trapezoid composed of 
6306 atoms (hgure 3). The curves for the 3NNTBA model have been horizontally shifted by 
AV^ = +0.072 Itnl/e as specihed in main text, section 2. Graph (a) show results for the electrodes 
attached at the 120° angles. Graph (b) shows analogous results for the 60° angles. The models of 
the electrodes are identical to those used for results in hgure 3. See also hgure 11 for a detailed 
view on the low-V^ interval of the present graph. 


(hgure 2) and the ZZ-terminated one composed of 6306 atoms (hgure 3). The comparisons are 
provided only for the AA-AA and OA-OA attachments. As opposed to main text, in this appendix 
we use the linear scales on the conductance axes. 
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Figure 11: A detailed view of the results from hgure 10. 
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Overall view of the results for the ac case is shown in hgure 8 where we display both the conduc¬ 
tance within the INNTBA (the solid green plots) as well as the results with the extended hoppings 
included (the dashed red plots). The conductance spectrum in the immediate neighborhood of the 
Fermi level is showed in a zoomed view in hgure 9. We see that the extended TB hopping range 
modihes the results only quantitatively. 

We now turn to the case of the ZZ-edge terminated trapezoidal hake. Results are shown in 
hgures 10 and 11. As we can see, in this case the extended TB hopping range has a noticeable ehect 
on the conductances within an about 0.02|fB|/e narrow gate voltage window. Most importantly, 
the AA-AA setup exhibits rapid oscillations of the conductance within this range [hgure 11(b), 
red plot], in this way opening the high-conductance window. Despite of this ehect, in most of the 
0.3 Itsl/e wide central gate voltage window the AA-AA setup remains insulating [hgure 10(b), red 
plot]. Similar rapid oscillations of the conductance are found also for the highly conductive OA-OA 
setup [hgure 11(a), red plot]. In this case the oscillations partially decrease the conductance within 
the narrow central window. Still, the conductance remains in average signihcant as it was found 
also within the INNTBA. 


B Electrodes composed of interacting wires 

As specihed in section 2, each of the two electrodes in our model is composed of a bunch of 
identical mutually non-interacting monoatomic wires; see hgure 1 which shows a schematic graphical 
representation of such an electrode and its contact to the graphene hake. It may be questioned 
whether such a model of the electrodes is sufficiently realistic. We note that in theoretical studies 
of electronic transport in GNRs very simple models of electrodes are often used, for example an 
electrode being just a continuation of the GNR [6, 7, 14, 17]. Despite of the substantial diherence 
between our model and the referenced ones we have shown [16], using AGNRs as examples, that our 
results are in semi-quantitative agreement with those employing the diherent models. We provide 
similar evidence also in SI, see section 1 and hgure 1 therein. 

In addition here we perform an explicit test how would our result change if the wires were 
mutually interacting. To achieve this we assume a model in which atoms of the wires within the 
20 a vertical distance from the hake surface become mutually coupled using the 3NNTBA model 
with the same parametrisation as used for the graphene hake. he. we assume an augmented central 
system formed by the hake and by the hnite pieces of the electrodes (FPE). The FPE are just those 
parts of the electrodes in which the wires are mutually coupled. All couplings within the augmented 
system are treated on equal footing using the 3NNTBA model. In this more complex model we still 
employ the semi-inhnite mutually decoupled wires. They are hxed to the ends of the FPE, not to 
the hake. The couplings within the FPE, although not describing any real system, provide us with 
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Figure 12: Linear conductance of the 6306-atoms large trapezoidal flake depicted on the left panel 
of hgure 3. The electrodes are attached at the 60° corners. Dashed black plot in the present hgure 
is just a detailed view of the black plot in hgure 3(b). Solid red plot in the present hgure describes 
the same hake contacted to the diherent model of the electrodes. Most importantly, the wires 
forming each of the two the electrodes are mutually coupled up to certain distance from the hake. 
See text in appendix B for more description. 


a verihcation example whether or not the internal structure of the electrodes is important for the 
ehects studied in the present paper. As the test case we choose the ZZ-edge terminated trapezoidal 
hake shown on the left panel of hgure 3, the second structure from the top. Results shown in 
hgure 12 clearly demonstrate that the diherent model of the electrodes has only a marginal impact 
on the results reported in our present work. We can make the conclusion that the ehects reported 
in the present paper can be found for various types of electrodes assuming that they have order of 
nanometer thickness. 


C Impact of the edge disorder 

In order to qualitatively assess the impact of the irregular edges on the main ehects reported in the 
present paper, we perform calculations for two representative junctions. First of them is based on 
the ac-edge terminated trapezoid shown in the left panel of hgure 2, the second structure from the 
top. After introduction of several random edge defect (including defects close to the corners) the 
structure consists of 6245 atoms (an image is provided in SI, hgure 7 therein). The perturbations are 
made in such a way that every carbon atom at the border has either two or three nearest-neighbour 
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Figure 13: Linear conductance of the sample with the edge defects. In the absence of the defects 
the flake would be identical to the one shown in the left panel of hgure 2, the second structure 
from the top. The electrodes are 34 monoatomic wires thick. Results for the AA-AA and OA-OA 
setups are shown only. Graphs with direct comparisons between the perfect and imperfect border 
cases and an image of the structure can be found in SI, hgure 7 therein. 

carbons. Similarly as for the structures with regular edges, here we also assume that the edges are 
terminated by hydrogen atoms [12]. All three bonding parameters tg and tg remain the same 
as for models with perfect edges. Similarly we keep all equilibrium on-site energies at zeros. Due 
to the defects present even at the corners the electrode’s thickness has been slightly reduced to 
34 monoatomically thin wires for each of the two electrodes. Calculated linear conductances for 
AA-AA and OA-OA setups are shown in hgure 13 and should be compared to hgure 2(b) which 
uses the same color coding. Graphs directly comparing the perfect and imperfect-border cases can 
again be found in SI, hgure 7 therein; linear scales on vertical axes are used there. The impact 
of the edge disorder is multiple: (i) the relatively uniform Gun prohle for the OA-OA setup now 
strongly varies with the gate voltage, (ii) the average magnitude of Gun has been decreased by 
about a factor of 10, i.e. the high conductance of the OA-OA setup has been lost, (hi) the Gun 
prohle for the AA-AA setup has got also a more oscillatory character, (iv) the Giin(I4) prohle for 
the AA-AA setup now exhibits signs corresponding to ZZ edges; see the central double peak in the 
black plot of hgure 13 just above the zero energy. The latter feature is not surprising because some 
of the defect in the ac edges provide pieces of ZZ-type termination. Despite of these modihcations 
the conductance of the AA-AA setup is in average still several orders of magnitude below those of 
the relatively conductive OA-OA setup. 

The second junction with imperfect edges is based on the ZZ-edge terminated trapezoid shown in 
the left panel of hgure 3, the second structure from the top. The hake with the defects is formed by 
6219 atoms, the electrodes are 33 wires thick and the other conditions and treatment are the same 
as in the ac case described above. The result are presented in hgure 14 and should be compared to 
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Figure 14: Linear conductance of the sample with the edge defects. In the absence of the defects 
the flake would be identical to the one shown in the left panel of hgure 3, the second structure 
from the top. The electrodes are 33 monoatomic wires thick. Results for the AA-AA and OA-OA 
setups are shown only. Graphs with direct comparisons between the perfect and imperfect border 
cases and an image of the structure can be found in SI, hgure 8 therein. 

hgure 3(b). As in the ac case, additional graphs as well as an image of the hake can be found in SI, 
hgure 8 therein. We observe that the disorder introduced to ZZ edges has a noticeable impact on the 
linear conductance function. This impact is smaller than in the ac case, fully in line with hndings 
of [12]. Most signihcantly, the narrow window of the ZZ-edge state induced high conductance (just 
at the Fermi level) becomes even narrower what is easily comprehended because the ZZ edges are 
now frequently interrupted by the defects. Away from the Fermi level the conductance of the OA- 
OA setup is also high and, interestingly, almost unahected by the disorder [compare blue plots in 
hgures 3(b) and 14]. At these gate voltages the Giin(14) profile of the (insulating) AA-AA setup is 
also little impacted by the ZZ edge defects [compare black plots in hgures 3(b) and 14]. 

We can hnally conclude that the main effects reported in the present paper are to a signihcant 
extent found operational also for hakes with edge disorder. 
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